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Abstract 

In order to better understand Kondo insulators, we have studied both 
the symmetric and asymmetric Anderson lattices at half-filling in one di- 
mension using the density matrix formulation of the numerical renormal- 
ization group. The asymmetric case is treated in the mixed valence regime. 
We have calculated the charge gap, spin gap and quasiparticle g ap as a 
function of the repulsive interaction U using open boundary conditions for 
lattices as large as 24 sites. We find that the charge gap is larger than the 
spin gap for all U for both the symmetric and asymmetric cases. RKKY 
interactions are evident in the f-spin-f-spin correlation functions at large 
U in the symmetric case, but are suppressed in the asymmetric the 
f-level approaches the Fermi energy. This suppression can also be seen in 
the staggered susceptibility x(q — 2/c/), and it is consistent with neutron 
scattering measurements of x(q) in CeNiSn. 

PACS nos.: 71.27.+a, 75.20.Hr, 75.30.Mb, 75.40.Mg 



1 Introduction 



Kondo insulators are a class of rare earth compounds that become semiconducting at 
low temperatures. However, they are not just ordinary semiconductors because their 
gaps are due to many body interactions between extended conduction and localized 
electrons in the f-shell. At high temperatures (T > 100K), these materials behave as 
metals with local magnetic moments, but at low temperatures they exhibit a small 
gap in their excitation spectrum and behave as narrow gap semiconductors (typically 
A ~ 100K). This class of materials includes CesBi4Pt3, CeNiSn, and SmBg [|TJ. It has 
also been argued that the transition metal compound, FeSi, is also a Kondo insulator 
i- 

Theoretically, the Kondo insulators have been modeled by both the Kondo lattice 
and the periodic Anderson model at half filling [fj, |3], f|, ||, |, 0. In the Anderson 
lattice each site has a localized orbital that hybridizes with an extended band of 
conduction electrons. Double occupation of the localized orbital is penalized by a 
strong Coulomb repulsion U. The half-filled system, in which there are two electrons 
per site, is an insulator at zero temperature. 

In the Anderson model, two distinct regimes can be considered. In the Kondo 
regime, the energy of the localized band lies well below the Fermi surface. In this case, 
the occupation of the localized level is very close to one, and charge fluctuations in and 
out of the localized orbital are negligible. On the other hand, in the mixed valence 
regime, the energy of the localized orbital lies very close to the Fermi level, and 
electrons can hop in and out of the localized orbital. As a result, charge fluctuations 
are important and the average occupancy of the localized orbital is less than one. 

Most of the theory that has been done on Kondo insulators has treated them in 
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the Kondo regime or assumed particle- hole symmetry or both |2|, ||, f|, |5], H Q]- (The 
symmetric Anderson lattice has particle-hole symmetry while the asymmetric case does 
not.) However, real materials do not necessarily have particle-hole symmetry, and it is 
generally believed that they are more likely to be in the mixed valence regime. Varma 
has argued that charge fluctuations are crucial to the understanding of the physics of 
these systems ||. He points out that such fluctuations should reduce the magnetic 
correlations between sites. As we shall see in section 4.2, this is confirmed by our 
calculations. Indeed neutron scattering experiments have indicated the absence of 
magnetic correlations in CeNiSn ||. In addition, measurements of the occupation of 
the f-orbital in Ce^Bi^Pt^ give rif = 0.865 at T = [fLOfl , indicating a mixed valence 
state. Thus it is important to determine the properties of the asymmetric Anderson 
lattice in the mixed valence regime. 

Previous studies of the mixed valence case have been hampered by various limita- 
tions. Analytic calculations, which include the Gutzwiller approach [TJ] and the 
slave boson mean field approximation [O, H|, do not distinguish between the different 
kinds of gaps (spin, charge and quasiparticle gap). Numerical approaches have been 



limited to 4-site chains [15, 16| and did not consider the charge gap 15 



In this paper we present the first systematic study of all three gaps of the half 
filled Anderson lattice in both the Kondo and the mixed valence regime. To the best 
of our knowledge, this is the first time that the staggered susceptibility x{q = 
the f-spin-f-spin correlation functions, and the occupation rif of the localized orbital 
have been calculated for the mixed valence case. 

The density matrix renormalization group technique allows us to consider chains of 
up to 24 sites. This gives us an advantage over other numerical approaches || [15], |TjJ 
which can only deal with short chains. Not only are we able to get f-spin-f-spin 



2 



correlation functions, but by considering long chains, we can also study the behavior 
of the f-spin-f-spin correlations in both regimes. We establish that the amplitude of 
the correlations decays exponentially in space, and we can calculate the correlation 
length as a function of the parameters. Previous Monte Carlo studies of the symmetric 
case were unable to distinguish between power law or exponential || because of the 
length of the lattice. In the mixed valence case, where the accuracy is higher than 
for the symmetric case at large U, we can get the extrapolated value of the spin gap 
and the average occupation of the localized orbital for infinite chains. We find that 
both have power law behavior for large U. This analysis was not possible in previous 



studies of short chains 15, 16 



We also consider the U = oo in which the doubly occupied states of the local- 
ized orbital are suppressed. Several approximation schemes, such as the Gutzwiller 
technique and the slave boson mean field approach, have been applied in this limit. 
However, it is not clear how well these approximations describe the properties of the 
Anderson Hamiltonian. To find out, we compare some of our results with those of 
the slave boson mean field approach. We find good agreement when the system is in 
the mixed valence regime but as the system approaches the Kondo regime, the slave 
boson picture breakdowns. 

The paper is organized as follows. We consider the one dimensional Anderson lat- 
tice at half filling using the density matrix formulation of the quantum renormalization 
group technique fli~7 , 18[] which is briefly described in section 2. After discussing the 
Anderson and Kondo lattice Hamiltonians in section 3, we give our results in section 
4. In both the symmetric and asymmetric case, we find in section 4.1 that the charge 
gap is larger than the spin gap. In section 4.2, we study the RKKY interactions 
by calculating the f-spin-f-spin correlation function and the staggered susceptibility 



3 



x(q = 2kf) as a function of U. Our results show that in the Kondo regime the RKKY 
interactions become important as U increases, while in the mixed valence case such 
magnetic correlations are suppressed because the f-orbitals are not always occupied. 
In section 4.3, we find that the occupation n/ of the localized orbital changes from 1 
at U = to a value close to 0.7 at U = oo in the asymmetric case when the localized 
band is at the Fermi energy. (In the symmetric case n/ = 1.) In section 4.4, we study 
the U = oo case. We present our conclusions in section 5. 



2 Quantum Renormalization Group Technique 

The density matrix formulation of the numerical renormalization group has proven 
extremely useful for the study of the ground state and low energy excitations of one- 



dimensional systems |fTP| , It has distinct advantages over other computational 
techniques. For example, it allows us to study chains longer than what is achievable 
by exact diagonalization, and it avoids the problem of Monte Carlo calculations where 
very low temperatures are difficult to attain. 

In real space renormalization group schemes, the system is divided into small 
blocks. The traditional renormalization group approach consists of diagonalizing the 
Hamiltonian of a small block, keeping the lowest energy eigenstates, and then forming 
a bigger block by combining a few of the small blocks. The drawback of this approach 
is that it completely neglects quantum fluctuations at the boundary of the blocks. 
For this reason, the standard approach performs poorly and there is not a significant 
improvement by increasing the number of states kept [[17 . 



The density matrix formulation of the renormalization group provides a systematic 
method of choosing the states to keep for the bigger block and at the same time it takes 
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into account the quantum fluctuations at the boundaries. The procedure for a one 
dimensional lattice goes as follows. We divide the chain into four blocks. Typically, 
blocks 1 and 4 have many sites and blocks 2 and 3 (in the center of the chain) are 
single sites. We find the ground state of the Hamiltonian for the lattice: 

>= ^ii,*a,<3,*4|*l, «2, «3, i4 > (1) 

where il, i2, i3, and z4 label states in blocks 1, 2, 3 and 4. Then we consider blocks 
1 and 2 to be the "system" and the rest of the chain to be the "environment". The 
objective is to construct a new block 1 out of blocks 1 and 2 by keeping only a fraction 
of the total number of states. In order to choose which states to keep, we form the 
density matrix for blocks 1 and 2 by tracing out the degrees of freedom of blocks 3 
and 4 [fZTJ: 



p(il, i2, zzl, U2) = J2 • ^tii,^,a (2) 

i3,i4 

We then diagonalize this density matrix. The eigenvalues of the density matrix give 
the weight of the associated eigenstates |^TJ. Since it is not possible to keep all the 
states, we truncate the basis by keeping a predetermined number of states with the 
largest weights. We use these eigenstates to construct the basis of the new block 1 
which is made out of the old blocks 1 and 2. In this way, we increase the size of block 
1 by one site. 

It is possible to target more than just the ground state; the lowest excited states 
can be calculated as well by simply diagonalizing the Hamiltonian for these states. In 
each case the diagonalization is performed in a subspace with fixed quantum numbers, 
e.g., with a fixed number of electrons and a fixed z-component of the total spin. 
This method performs at its best for open boundary conditions in which there is no 
hopping past the ends of the chain. We typically keep 100 states, although for cases 



where higher accuracy is needed, we keep up to 200 states. For small values of the 
Coulomb repulsion, the results are very accurate with truncation errors of the order 
of 10 -8 , while for large values of the Coulomb repulsion the accuracy is reduced and 
the truncation errors increase to 10~ 4 . For a more detailed description of the method, 
see reference [17]. 

As a check for our method, we compare our energies and correlation functions 
with exact diagonalization results. For short chains we can keep all the states, and 
we obtain full agreement with exact diagonalization. As another check, we compare 



the ground state energy for the symmetric Anderson lattice with Gutzwiller [2^] and 
Monte Carlo results @ as shown in Fig. 1. We plot (e D + Z7/2)/|e |, where e Q is the 
ground state energy per site. The Monte Carlo simulations calculated the ground 
state energy for 16 site chains with periodic boundary conditions jTJ. We use open 
boundary conditions and chains of 24 sites. The Gutzwiller values are for a correlated 



spin density wave ||22|| . We obtain good agreement with both calculations. For large U, 
our energies lie below the Monte Carlo and the Gutzwiller values, giving a better upper 
bound for the ground state energy. In this regime the binding energy is proportional 
to l/U in agreement with perturbation theory. 



3 The Periodic Anderson Hamiltonian 

We consider the standard periodic Anderson Hamiltonian in one dimension: 

H = ~t J2( C i*Ci+lv + 4+1*0 + 6f J2 n L + U J2 n Wil + Vj2( C laU + flcia) (3) 
icr itr i icr 

where c\ a and Cj CT create and annihilate conduction electrons with spin a at lattice site 
i, f\ a and fi a create and annihilate local f-electrons, t is the hopping matrix element 
for conduction electrons between neighboring sites, e/ is the energy of the localized 
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f-orbital, U is the on-site Coulomb repulsion of the f-electrons, and V is the on-site 
hybridization matrix element between electrons in the f-orbitals and the conduction 
band. For simplicity we neglect orbital degeneracy. U, V, t, and e/ are taken to be 
real numbers. 

Let us examine the Anderson Hamiltonian in various limits. For U — this 
Hamiltonian can be exactly diagonalized in k-space. We obtain two hybridized bands 
with energies : 

1 



± 



A, 



(e f - 2tcos(ka)) ± J(e f + 2tcos(ka)) 2 + AV 2 



(4) 



where a is the lattice constant. When there are two electrons per unit cell, the lower 
band is full while the upper one is empty. Thus the system is insulating when N e = 2L. 
Here N e is the number of electrons and L is the number of sites in the lattice. 

When the mixing term is small (irV 2 /2t(ef + [/) « 1 and irV 2 /2tef < 1), the 
Anderson Hamiltonian can be mapped into the Kondo model |23| 



H K = -tJ2(4aCi+la + c\ +1<J C ia ) + J ef f ^ S{ • ^ (5) 
ia i 

where is the spin density of the conduction electrons at site i and J e ff is given by 



the Schrieffer- Wolff transformation 23 



2\V\'U 

Note that for the symmetric case where e/ = —U/2, J e ff = 8V 2 /U, and for U = oo, 
Jeff = ~2V 2 /e f . 

Although there is no direct interaction between the f-electrons in the Hamiltonian 
( 0), there is an indirect exchange coupling between the f-sites via the conduction band. 
This is the RKKY coupling, originally derived for the magnetic interaction between 



nuclei in metals [El]. An electron in the f-orbital tends to polarize the spin of the 



conduction electron cloud around it. The spin of the conduction electron polarization 
cloud oscillates with distance and it polarizes the spin of the f-electrons at neighboring 
sites. In the case of the Anderson Hamiltonian, the RKKY coupling constant can be 
obtained by considering two Anderson impurities in a gas of conduction electrons. The 
hybridization term is taken as a perturbation, and the fourth-order contribution to the 
perturbation expansion gives the coupling between sites. The effective Hamiltonian 
for the two impurities is: 

Hrkky = JrkkyS[ ■ S 2 (7) 

where S( and £2 are the spin operators for the f-electrons on the two impurities. The 
explicit expression for Jrkky is fairly complicated in general , but it can be shown 
that it favors antiferro magnetic alignment of the f-spins with wavevector q = 2k f, 
where kf is the Fermi wavevector of the noninteracting conduction electrons. 

For the one dimensional symmetric case (ej = —U/2) at half filling, 2k^a = n, and 
the asymptotic expression for large distances is |J 

jRKKY -{ir) 8^~ (8) 

— * — * 

where / is the distance between Si and S2 in units of the lattice constant. Note that 
Jrkky = (— l)'Jf/y/87rf/, that is, the Kondo coupling constant is the parameter that 
governs the f-f correlations in the Kondo regime. The asymptotic expression ( [|) is 
the same as the one obtained by starting with the Kondo Hamiltonian and doing 
second-order perturbation theory in J e ff- 

We expect RKKY interactions to be important in the Kondo regime where the oc- 
cupancy of the f-orbital is very close to 1. In the mixed valence regime, the occupancy 
of the f-level is reduced, and therefore RKKY correlations should be suppressed |§. 
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3.1 Symmetric case 

The symmetric case corresponds to a value of 6f = —U/2. For large U, the f-level 
lies well below the Fermi surface and the system is in the Kondo regime. Particle- 
hole symmetry ensures that the occupancy of the f-level, nf, is always one. In the 
strong coupling regime (8V 2 /U <C t), charge fluctuations in the f-orbital are highly 
suppressed because the Coulomb repulsion makes the states with one electron in the 
f-level much more energetically favorable than the state with two electrons or the state 
with the orbital empty. 

In addition to the usual SU(2) spin symmetry, this particular choice of e/ intro- 
duces an SU(2) charge pseudospin symmetry into the system ||. The pseudospin 
operator (/) can be obtained from the spin operator S by performing a particle-hole 
transformation on one of the spin species ||. Its components are: 

h = lE( c h^+4 l c ii +fy i1 +fif il -2) 

i 

u = E(-!) 1 '(4t4i - fyio (9) 

i 

i- = E(- 1 )'( c ^ c n - hfo) 

i 

The z component of the pseudospin is equal to N e i/2 — L. Note that half filling 
corresponds to I z = 0. An I z = 1 state can be achieved by adding two electrons. 

All the energy eigenstates of the symmetric Anderson model have a definite value 
of S and I. At half filling the ground state is a singlet both in spin and pseudospin 
space (S = 0, / = 0) PS| . The lowest lying excited state is a spin triplet. The energy 
difference between these two states is the spin gap A s : 

A s = E{S = 1,1 = 0) - E (S = 0,1 = 0) (10) 

where E Q is the energy of the ground state. 
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To find the charge gap, we note that optical experiments measure the charge gap 
by measuring the conductivity which is determined by the current-current correlation 
function. The current is related to the charge density through the continuity equation. 
Thus the lowest lying charge excitation is the lowest excited state \n > with S = 
such that < 0| J2 q Pq\ n >7^ 0, where p q is the q-component of the Fourier transformed 



charge density operator and |0 > is the ground state p6| . Notice that p q is related to 



Ig, where I q is a Fourier transformed vector in pseudospin space given by 

i 

j- = ^-^(-1)%^ - f n f n ) 

i 

Using the Wigner-Eckart theorem, one can show that the (1=1, S=0) states are the 
only states \n > for which the charge density p q has finite matrix elements < n|p 9 |0 > 
with the ground state |0 >. Thus the charge gap A c is the energy difference between 
the ground state and the lowest pseudospin triplet state |J: 

A c = E(S = 0,1 = 1) - E (S = 0,1 = 0) (12) 

The quasiparticle gap gives the energy for making a noninteracting particle and 
hole. It is defined as follows: 

A qp = E (N e = 2L + 1) + E (N e = 2L - 1) - 2E (N e = 2L) (13) 

Due to the particle-hole symmetry of the symmetric case, one can write 

A qp = 2[E (N e = 2L + 1) - E {N e = 2L)\ (14) 

The quasiparticle gap can be thought of as the difference of chemical potentials 

Agp = fI N+1 - p N (15) 
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where /i N = E Q (N e ) - E Q (N e - 1) and N e = 2L for half filling. 

For U = 0, all the gaps coincide and are given by the hybridization gap: 

A s = A c = A qp = \+-\- = 2Vt 2 + V 2 - 2t (16) 

As U — > oo, the f-electrons decouple from the conduction electrons, and all the gaps 
go to zero. 

We can also define the gap A ns between the ground state and the lowest excited 
neutral singlet state which has quantum numbers (S=0, 1=0) 0. For the half filled 
Kondo lattice, the neutral singlet has been found to be an elementary excitation 
consisting of a "particle" and a "hole" which are (S=l/2, 1=1/2) excitations. In a 
single site basis, a "hole" is a site with one f-electron and no conduction electrons 
with quantum numbers (S=l/2, 1=1/2, l z = —1/2), while a "particle" is a site with 
one f-electron and two conduction electrons with (S=l/2, 1=1/2, l z = +1/2). The 
Kondo and Anderson lattices agree for small J e ff/t. In this regime the neutral singlet 
gap is smaller than the charge and quasiparticle gaps but larger than the spin gap 0. 
However, for small U, the lowest lying (S=0, 1=0) excited state is not an elementary 
excitation, but presumably consists of two spin excitations. Indeed, for U — 0, one 
can show that A ns = 2A S . In the small U limit, there are bands of spin and charge 
excitations which lie between the lowest excited neutral singlet state and the ground 
state. Because of this, it is hard for the numerical renormalization group technique 
to accurately calculate the energy of the neutral singlet state for small U. For this 
reason, we have not calculated the neutral singlet gap. 
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3.2 Asymmetric case 

The symmetric case is typically studied in the strong coupling regime (8V 2 /U <C t) 
where the f-level lies far below the Fermi energy. However, in real materials the f-level 
can be near the Fermi surface even though U is large. In addition real materials ordi- 
narily do not have particle-hole symmetry. For this reason we consider the asymmetric 
case in which e/ can have any value irrespective of U. In particular, we place the f- 
orbital right at the Fermi level (e/ = 0) so that states with no electrons and states 
with one electron in the f-level are equal in energy, while states with two electrons are 
highly suppressed due to the strong Coulomb repulsion. In this mixed valence case, 
charge fluctuations are allowed and the occupation of the f-level, rif, is always less 
than one for U > 0. 

Note that there is no particle-hole symmetry and the pseudospin operator is no 
longer conserved. Only the total number of electrons (I z ) remains a good quantum 
number. For all the values of the parameters that we explored, we found that the 
ground state is still a spin singlet at half filling. There is still a gap to the lowest lying 
excited state which is a spin triplet. Thus the spin gap can still be defined as: 

A s = E(S = 1) - E (S = 0) (17) 

The quasiparticle gap is defined in equation ( [13]). Note that in this case E(2L + 1) ^ 
E(2L — 1) since there is no particle hole symmetry. 

Unlike the symmetric case we can no longer use charge pseudospin symmetry to 
find the charge gap. We must use the fact that the lowest lying charge excitation is 
the lowest excited state |n > with S = such that < n\ Z) 9 p 9 |0 0, where p q is 
the q-component of the Fourier transformed charge density operator and |0 > is the 
ground state [^j. Using the fact that J2 q Pq — c !=o c «=(b we find that |n > is the lowest 
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S = excited state. Let its energy be denoted by Ei(S = 
given by 

A c = E 1 (S = 0)-E o 



0). Then the charge gap is 



(18) 



4 Results 

In using the density matrix renormalization group technique, we consider lattices up 
to 24 sites long, and we use open boundary conditions. We set t=l so that energies 
are measured in units of t. For the symmetric case, by definition e/ = —U/2. In order 
to study the asymmetric case in the mixed valence regime, we set 6f = 0, that is, tf is 
at the Fermi energy. We fix V=l and vary U from to 30. Although this value of V 
is somewhat larger than in real materials, it is convenient for numerical reasons and 
for comparison with previous calculations ||. 

4.1 Gaps 

In Figs. 2a and 2b, we plot the gaps versus U for different chain sizes for the symmetric 
and asymmetric cases, respectively. 

Relative Gap Sizes. In both the symmetric and asymmetric cases, the charge and 
quasiparticle gaps are larger than the spin gap for U > 0. This is consistent with 
experiments which find that the charge gap is larger than the spin gap in Ce3Bi 4 Pt 3 
f2"8fl . As we shall see, in the mixed valence case the ratio A c /A s ~ 2 for large U. This 
agrees with the experimental value of 1.8 for CesBi4Pt3 |[28|| . In the Kondo regime, 
one expects a much larger ratio when U is large; in fact, this ratio diverges in the 
symmetric U — > oo HI. 

Note that the gaps have a smaller value at U = oo than at U = 0. In the 
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asymmetric case this reduction is not as large as in the symmetric case. This decrease is 
due to many body interaction effects and may explain why band structure calculations 
p9j |30| find a larger gap than the optical gap measured experimentally |31 . 



Nature of the Excitations. For the symmetric case (Fig. 2a), notice that the spin 
gap decreases monotonically with increasing U. In contrast to the monotonic behavior 
of the spin gap, Fig. 2a shows that the charge and quasiparticle gaps initially increase 
with U, go through a maximum and then decrease. This behavior was also observed 



by Steiner et al. f2"T|| . To understand this, note that the quasiparticle gap is obtained 
by adding a particle to the system. For very small U, this particle goes predominantly 
into an f-orbital; so as U increases, the gap also increases due to the Coulomb repulsion. 
However, for large U, the extra particle goes mostly into the conduction band, and 
the gap starts decreasing. One can see how much of an additional particle goes into 
an f-orbital by plotting < Nf > —L vs. U for the state with N e i = 2L + 1 (see Fig. 
3). Since the total number of f-electrons, Nf, equals L at half-filling with particle-hole 
symmetry, < Nf > —L tells us how much of the extra particle is in the f-orbital. One 
can see in Fig. 3 that the extra particle goes mainly into an f-level for small U but, 
as U increases, it goes more and more into the conduction band. 

Now we discuss the charge gap shown in Fig. 2a. In the symmetric case, the 
charge gap is obtained by adding two particles to the half-filled system to get an 
(J = 1, I z = 1) state. It has a maximum as a function of U for the same reason that 
the quasiparticle gap does. The fact that A c > A qp for any finite value of U means 
that the two extra particles repel each other. So as L — > oo, they will be infinitely 
apart and A c will equal A qp . 

In Fig. 2b we show the gaps versus U in the asymmetric case. We see that for 
U ~ 2, A c is greater than A qp , but as U increases, they cross. This crossover can be 
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interpreted as follows. The quasiparticle gap consists of two noninteracting 5=1/2 
excitations. The charge gap, on the other hand, is given by the lowest S = excited 
state. For small U this state consists of two 5 = 1/2 excitations, while for large U 
there is a crossover to a state made out of two S = 1 excitations. This picture is 
supported by the fact that for small U, A c ~ A qp , but as U increases, A c ~ 2A S . 

Finite Size Effects as U — > oo. As U — > oo in the symmetric case, the conduction 
electrons and the f-electrons decouple, and the gaps approach their values for a free 
electron band on a finite size lattice with open boundary conditions. We confirmed 
this by calculating the gaps on finite size lattices for free electrons. The charge and 
quasiparticle gaps have large finite size effects, especially for large U. Nishino and 
Ueda H find a much smaller charge gap in the large U limit because they use periodic 
boundary conditions, and the charge gap goes to zero for free electrons on a ring. On 
the other hand, our values for the spin gap agree with the results of Nishino and Ueda 
H because the spin gap has small finite size effects. To understand this, note that 
the the spin excitation has small dispersion because it is a spin flip that has mainly 
f-character. This makes it a local excitation which is not very sensitive to the length 
of the lattice or to the boundary conditions. To confirm this picture, in Fig. 3 we plot 
< S£ > versus U for the excited (S — 1, 1 = 0, S z — 1) state. (Here S£ refers to the 
z-component of the total f-spin of the chain.) For U = 0, < S{ > is about 80% of the 
total spin and it gets very close to 100% for larger U. 



In the asymmetric case, the gaps tend to a finite value |L5] as U — > oo. The fact 
that the gaps are finite is not a finite size effect, but rather can be understood in the 
following way. For U = oo, states with double f-occupancy are suppressed, but the 
proximity of the f-level to the Fermi energy allows the system to fluctuate between 
states with no f-electrons in the local orbital and states with one f-electron. This 
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leads to hybridization and the formation of the gap. Another way to see that there is 
hybridization is to recall from eq. ( |6|) that J e ff = —2V 2 /ef for U = oo. The presence 
of a finite J e // leads to hybridization and gap formation. 

Spin Gap. For the asymmetric case, we can study how the spin gap approaches its 
U = oo value. We begin by extrapolating the spin gap to its value at L — > oo in the 
following way: for each value of U, we plot A s versus 1/L 2 , and obtain a straight line 
(Fig. 4). The intercept is then the value for L — > oo. To understand why the spin gap 
is proportional to 1/L 2 , note that the spin density versus site for the (S = 1, S z = 1) 
state is shaped like a particle-in-a-box wave function of wavelength ~ L. This is shown 
in the inset of Fig. 4. We then interpret the spin gap as the kinetic energy of the 
quasiparticle which is proportional to the square of the wavevector and goes as 1/L 2 . 

We study the behavior of the spin gap for L — > oo and we find that, for large U 
with ef = in the asymmetric case, it can be fitted by the form 

A S (U) = A S (U = oo) + (19) 

with a — > 1 as U — > oo. For e/ < 0, a study of the large U behavior of the gap for 
short chains suggests that the spin gap approaches a finite value as (ef + U)^ 1 for the 
asymmetric case. 

4.2 RKKY Interactions 

We study RKKY interactions between f-spins in the ground state by calculating the 
f-spin-f-spin correlation function and the staggered susceptibility. We first turn our 
attention to the correlation function. 
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4.2.1 f-spin-f-spin Correlations 

We study the f-spin-f-spin correlations in the ground state due to RKKY interactions. 
In Fig. 5 we show the spin-spin correlation function of the f-electrons as a function of 
spatial separation r (measured in units of the lattice constant) for different values of 
U. Fig. 5a corresponds to the symmetric case and Fig. 5b to the asymmetric case. 

Mixed Valence Case. In the mixed valence case, the RKKY correlations are highly 
suppressed because the f-orbital is not always occupied. This can be seen in the f- 
spin-f-spin correlation function shown in Fig. 5b. In contrast to the symmetric case, 
we see that the increase of U by an order of magnitude does not change the rapid 
decrease of the correlations with distance. 

It is tempting to speculate that another reason why the RKKY correlations are 
suppressed is that the effective Kondo coupling J e ff becomes very large as ej ap- 
proaches the Fermi level. (Recall that J e ff = —2V 2 /tf for U = oo.) In the Kondo 
lattice model, RKKY correlations decrease with increasing J e ff- The problem with 
this explanation is that the Schrieffer- Wolff transformation is no longer valid when 
|e/| <C T where T = 7cV 2 p. Even with more sophisticated techniques |33], there is no 
known analytic expression for J e ff when \ef\ <C Y and Y/U <C 1. 

Symmetric Case. For the symmetric case we see that as we increase U from 2 
to 20, the amplitude of the oscillations becomes much larger and persists for longer 
distances. There are several ways to understand why RKKY interactions increase with 
large U. First one notes that as U increases, the effective Kondo coupling constant 
J e ff = 8V 2 /U decreases relative to the hopping matrix element t. Increased hopping 
means increased correlations between sites. The decrease of the gaps also enhances 
RKKY interactions. Another way to understand the growth of RKKY interactions is 
to note that the Kondo temperature T% ~ exp(— l/pJ e //) decreases with increasing U 
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and, as a result, the Kondo compensation cloud increases in size. Uncompensated f- 
spins within the cloud develop correlations || . A third way to view this was suggested 



by Varma and Doniach [fffil . In the Kondo model they have argued that RKKY 
interactions will dominate as J e ff decreases because the RKKY energy scale goes 
as Jgjj while the Kondo energy scale depends exponentially on J e ff- Opposing this 
enhancement of RKKY interactions is the fact that the effective RKKY coupling will 
decrease as J e ff decreases, but only as Jiff- 

We find that the amplitude of the RKKY oscillations decreases exponentially with 
distance, and it can be fitted by exp(— r/£), where £ is the correlation length. In Fig. 
6, we show the behavior of the correlation length for the symmetric function 
of the effective Kondo coupling constant J e ff = 8V 2 /U. We plot two different sets of 
parameters and we see that the curves lie on top of each other, confirming the fact 
that the relevant parameter in the strong coupling regime is J e ff- We also show the 
results for the one dimensional Kondo chain, which agree very well for small J e ff- As 
J e ff decreases (i.e., as U increases), the correlation length goes up to a maximum and 
then starts going down. There are two ways to interpret this. One view says that 
this is to be expected since the f-electrons decouple from the conduction electrons as 
J e ff goes to zero, and the correlations between different sites decreases. On the other 
hand, one can interpret the exponential decay of the spin correlations with distance 
as a consequence of the existence of a gap in the excitation spectrum. The gap 
suppresses low energy (long wavelength) magnetic excitations, and this greatly reduces 
long range correlations. As the gap decreases, one would expect the correlation length 



to monotonically increase roughly as £ ~ vp erm i/ A$ [|l4 j. Since we find a maximum in 
£, our data supports the first point of view, although the error bars on £ are bigger for 
small J e ff because in this regime the effects of the ends affect the wavefunction of the 



18 



entire lattice. This can be seen, for example, in the oscillations of the hybridization 
energy from site to site. When J e ff is small, these oscillations remain sizable even in 
the center of the lattice. For large J e //, these oscillations decrease rapidly as one goes 
away from the ends of the chain. 

4.2.2 Staggered susceptibility 

In order to further study the magnetic properties of the system, we calculate the 
staggered susceptibility x(q) as a function of the momentum q. To calculate x(q)i we 
apply a small magnetic field h = h Q cos(gr) in the z direction for different values of 
h and plot S(q) versus h . S(q) is the Fourier transform of < 0|S' z (i)|0 >. (The field 
couples to both the f-spins and the conduction spins.) If h is small enough, this plot 
is a straight line whose slope is the susceptibility. For this application, q has to be a 
good quantum number. Therefore we use periodic boundary conditions. This greatly 
reduces the accuracy of the method, and we can only study short chains (4 or 6 sites). 

In Fig. 7 we plot the staggered susceptibility x(<? a — 2k fa = ir) versus U . We 
see that in the symmetric case the staggered susceptibility diverges as U — > oo. This 
is consistent with the growth of RKKY oscillations. In addition, as U — > oo, the 
f-electrons decouple from the conduction electrons, and become polarized in an arbi- 
trarily small magnetic field. As a result, the susceptibility will diverge for all q. In 
contrast, in the mixed valence U — > oo, the staggered susceptibility x(q — 2&/) 

approaches a finite value x(U = oo). x(U) does not change much with U because of 
the reduction of the magnetic correlations in the mixed valence regime. 

In the inset of Fig. 7, we show x{o) versus q for both the symmetric and asymmetric 
cases. In the asymmetric case the peak at q = 2k f is greatly reduced, that is, x(q) 
shows little structure. This is consistent with neutron scattering results on CeNiSn 
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U in which they find that x(q) * s independent of q within experimental error. In fact 
the fluctuations we find in x{q) are an order of magnitude smaller than those seen 
experimentally, though one should keep in mind that our lattices are small and the 
values we chose for the parameters are not necessarily appropriate for CeNiSn. 

In the mixed valence regime the occupation of the f-level is less than 1 per site, 
which implies that the conduction electron occupation is greater than 1 per site. One 
might worry that this implies that the RKKY wavevector 2k f is no longer equal to tt, 
but rather is equal to a wavevector that is incommensurate with the lattice. Could 
this explain the lack of divergence for x{<l a = 7r ) seen m Fig. 7? The answer is no. The 
suppression of RKKY correlations can clearly be seen in Fig. 5b. Since Fig. 5b is in 
real space, any RKKY correlations, even those with an incommensurate wavevector, 
would show up in this plot. 

4.3 Occupancy of the f-level 

In the symmetric case, the occupation of the f-level is 1 in the ground state due 
to particle-hole symmetry. However, in the asymmetric case we expect the f-level 
occupation rif to be less than 1 in the mixed valence regime. We can calculate the 
f-level occupation rif per site by taking the expectation value of nji on each site with 
respect to the ground state wavefunction and then averaging over the sites. Thus rif 
is given by 

™/ = 7E<°l/U<x|0> (20) 

L ia 

where |0 > is the ground state. In Fig. 8, we plot rif(L = oo) versus U. We can 
extrapolate to infinite lengths because rif versus 1/L follows a straight line due to 
the factor of 1/L in eq. ( |2D[) We see that as U increases, rif decreases from 1 
towards a value close to 0.7, indicating a mixed valence state at large U. For large U, 
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rif follows power law behavior: 

n f (U)= n f (U = oo) + (21) 

where n f (U = oo) = 0.675 ± 0.001 and (3 goes to 1.4 ± 0.05 as U -»• oo. 

The results shown in Fig. 8 are for e/ = 0. We have studied short chains to 
see what happens as drops below zero. In the mixed valence regime, our study 
suggests that rif has the same form as ( |21]) for large U. As e/ becomes more negative, 
(3 increases and approaches a value close to 2 for large \ef\. We show in Appendix A 
that perturbation theory agrees with this result for V 2 /\ef \ <C t. 



4.4 U = oo 

Our renormalization group approach allows us to study the U = oo case by eliminating 
the doubly occupied f-states from the Hilbert space. This allows us to compare our 
results to the predictions of approximation schemes which are applied in this limit, 
e.g., the slave boson approach |13| and the Gutzwiller approximation |TT |. 



We study what happens as we reduce V in the cases when e/ = and e/ = —0.5. 
In Fig. 9 we compare the spin gap with the one predicted by the slave boson mean 
field technique |Tj|. For the case e/ = shown in Fig. 9a, there is overall qualitative 
agreement, although the slave boson mean field theory overestimates the size of the 
gap. It is important to point out that the slave boson curve has been calculated using 
a constant density of states while the numerical results use a tight binding density 
of states. Nevertheless, when e/ = 0, both give power law behavior for small V | |15| |: 



A s ~ V A . For the case Ef — —0.5, the situation is quite different as can be seen in Fig. 
9b. There is qualitative agreement for V greater than 0.5 but then the slave boson 
curve drops very rapidly for small V. A similar situation can be seen in Figs. 10a and 
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10b for 1 — rif. 

To understand these results let us take a closer look at the slave boson mean field 
solution [[R|. The number of slave bosons J2ib\bi represents the number of sites with 



no f-electrons. Each site is subject to the constraint 

blk + £ ILL,, = 1 (22) 

Since the slave boson number operator is positive definite, each site is constrained to 
have less than one f-electron. In the mean field approximation the boson operator 
is approximated by its mean value a , i.e., b\ — — a Q , where a is a real number. 
This amounts to suppressing double occupancy of the f-level only on average. Using 
slave bosons, an effective Hamiltonian can be written in which doubly occupied sites 
have been projected out. As a result, there is no U-term. The resulting effective 
Hamiltonian is then equivalent to the U = Anderson Hamiltonian ( |3|) with an 
effective hybridization 

V? ff = V 2 al = V 2 (l-n f ) (23) 

The slave boson technique recovers the expected picture of two hybridized bands with 
a gap. The expression for the gap is the same as for the U = Anderson Hamiltonian 
with this effective hybridization. For small V, the gap goes as V^Jt. The f-level 
is also renormalized to if, and self-consistent equations for both if and a Q can be 
obtained [T3[. The small V expansion of the self-consistent equations gives 



£ t = 2ta 2 



o 



When 6f < and V is small enough such that a Q is very small, then if in the 
exponential can be neglected, and a D goes to zero as exp(— 2t\ef\/V 2 ). (if > 0.) 
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Thus the gap will depend exponentially on V. When e} becomes greater than e/, 
a 2 = 1 — rif crosses over from exponential to power law behavior. If |e/| <C e}, then 
e/ in the exponent of equation ( |24|) can be neglected. Plugging in e} = 2ta 2 Q gives a 
transcendental equation for a a . Solving this numerically, we find that a 2 ~ V 2 . As a 
consequence, the gap goes as V A . Thus the slave boson mean field approach predicts 
that the gap crosses over from exponential dependence on V to power law dependence 
as e/ approaches from below. However, our numerical renormalization group results 
do not show such a change in behavior. Although it might be argued that finite size 
systems will give power law rather than exponential behavior, increasing the length 
of the chain does not indicate any evidence of such a dramatic crossover. In addition, 
perturbation theory for an infinite chain shows that 1— nf goes as V 2 when V 2 /ef <C t 
(see Appendix A). This is exactly what we find for tf = —0.5 with small V. 

The reason for this disagreement between the slave boson calculation and the 
numerical renormalization group approach may be that the mean field approximation 
neglects spin fluctuations. Spin fluctuations are crucial in the Kondo regime. They 
also give rise to RKKY interactions. 



The U — oo limit has also been treated by the Gutzwiller approximation |11 



It predicts a ferromagnetic ground state when the system goes towards the Kondo 
regime where V 2 /tft <C 1. However, Rice and Ueda [I I] only considered uniform 
magnetic states, as opposed to antiferromagnetic states which are favored by RKKY 
interactions. Using a saddle point formulation of the Gutzwiller method, Reynolds 



et al. |35| suggest that the Gutzwiller approximation is biased too much in favor 
of a ferromagnetic state. However, they still find a ferromagnetic instability in the 
Kondo regime. We do not see any evidence of ferromagnetism in any of the cases we 
studied. We always find an S = ground state. We find that strong antiferromagnetic 
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correlations appear as the system goes towards the Kondo regime. 

5 Conclusions 

We have studied the one dimensional Anderson lattice at half filling using the density 
matrix numerical renormalization group technique. We considered the symmetric 
Anderson model with 6f — —U /2, as well as the asymmetric case in which we set = 
in order to study the mixed valence regime. Since many real materials lack particle- 
hole symmetry and are likely to be in the mixed valence regime, we concentrated on 
the asymmetric model with tf = 0. Our results suggest that some Kondo insulators 
behave as mixed valence systems, and therefore should be modeled by the Anderson 
lattice model in the mixed valence regime as opposed to the Kondo lattice or the 
Anderson lattice in the Kondo regime. 

For all the values of the parameters that we have used, including ej < in the 
asymmetric case, we find that this is an insulating system with gaps. As U — > oo in the 
symmetric case, the conduction electrons decouple from the f-electrons and the gaps 
approach the values that they have for free electrons. However, in the asymmetric case 
the gaps approach a finite value in the large U limit because the conduction electrons 
do not decouple from the f-electrons. There is hybridization due to charge fluctuations 
in and out of the band of local orbitals whose energy is at or near the Fermi energy. 

In both the symmetric and asymmetric cases we find that the charge and quasipar- 
ticle gaps are larger than the spin gap for U > 0. However, we found that the relative 
values are different in each case. For the mixed valence case the ratio between the 
charge gap and the spin gap is roughly two in the strong coupling regime (large U). 
This is in agreement with reference [EHJ in which the authors report a ratio between 
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the charge gap and the spin gap of Ce^Bi^Pt^ equal to 1.8. In contrast, the strong 
coupling limit of the symmetric case (Kondo regime) gives a much larger ratio between 
the two gaps. In fact, it was shown by Nishino and Ueda || that this ratio diverges 
as the Coulomb interaction U goes to infinity. Our result, together with the fact that 



the occupation of the localized orbital is 0.865 at T=0 [10], suggests that Ce^Bi^Pt-^, 
behaves like a mixed valence compound at low temperatures. 

In the asymmetric case the nature of the lowest lying charge excitation changes 
as U increases. For small U the charge state is made out of two spin-1/2 excitations, 
while for larger U, it becomes a state consisting of two spin-1 excitations. 

The fact that the gaps are smaller for large U than for U = is consistent with 
the fact that band structure calculations on FeSi [ESL BQ], which ignore many body 



correlations, predict values for the gap that are larger than the optical gap measured 
experimentally [[3l[ . 



We have studied the RKKY interactions by calculating the f-spin-f-spin correlation 
function as well as the staggered susceptibility x(l)- We find that the amplitude of 
the RKKY oscillations decays exponentially with distance. For the symmetric Ander- 
son model, the RKKY interactions become important in the strong coupling regime 
(8V 2 /U <C t). This can be seen, for example, in the divergence of the susceptibility 
x(q = 2kf) as U — > oo. This stands in sharp contrast to the asymmetric case in the 
mixed valence regime where the RKKY correlations are strongly suppressed due to the 
reduction in the f-occupancy. This agrees with Varma's argument || that magnetic 
correlations are suppressed in the mixed valence case. Our findings are also consis- 
tent with neutron scattering results on CeNiSn || which found x(o) independent of 
q. Further experiments would be necessary to confirm that CeNiSn is a mixed valence 
compound. One such measurement would be photoemission which can measure the 
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occupation of the f-level. In the mixed valence case that we consider (the localized 
level right at the Fermi energy and T ~ 1), we find that rif varies from 1 at U = to 
a value close to 0.7 at U — oo. 

We would like to thank Steve White, Herve Carruzzo, Andy Millis, Luiz Oliveira, 
Doug Scalapino, George Gruner, Kazuo Ueda, and John Wilkins for helpful discus- 
sions. This work was supported in part by ONR Grant No. N000014-91-J-1502 and 
an allocation of computer time from the University of California, Irvine. C. C. Yu is 
an Alfred P. Sloan Research Fellow. 

5.1 Appendix A 

In this appendix we use perturbation theory to show that for small V, 1 — rt/ ~ V 2 
for U = oo. We also show that nf(U) — Uf(U = oo) ~ U~ 2 . Consider the Anderson 
Hamiltonian ( ||) for the case where the hybridization matrix element is zero, i.e., 
V — 0, ef < 0, and €f + U > 0. In this case the ground state consists of a half-filled 
system with one electron per site in the f-band. The f-electrons states are degenerate 
with respect to their spin configurations. 

For V > but V 2 /2t(e f + U) < 1 and V 2 /2te f < 1 , perturbation theory can 
be used to calculate the correction to the ground state energy. The first non-zero 
contribution is the second order term: 

I < m\H AO > I 2 . . 

A£ = E ' J_ L ' (25) 

where H\ is the hybridization term. The only possible states that contribute to the 
sum are those that have one hole in the Fermi sea and an extra electron in the f-band 
or vice- versa. Also, the fact that kfa = ir /2 = ir — k^a allows one to write these two 



26 



contributions under the same sum: 



A£ = V 2 Y, 

k>kf 



1 



1 



(26) 



_e/-e fe e k + e f + U_ 

where e k is the energy of the conduction band. In our case, e k = —2tcos(ka). For 
the symmetric case, e/ = —U/2, this expression coincides with the one obtained by 
Blankenbecler et al. J7|. The total ground state energy to second order is given by 

E n 



, Y,k<k f e k AE Q 



AT 



(27) 



The occupation of the f-orbital can be obtained by the relation 



For U — oo, this expression gives: 



+ 



1-71,(^=00) = — £ 



_(e fc + e/ + [/) 2 (e/-e fc ) 



I' 2 



k>k f v fc / 



( e / - e fc ) 



This implies that for small V, 1 — rif ~ V 2 . Note that for arbitrary U 



n f{U) - n f (U = oo) 



N 



E 



k>k 



' (e k + e f + Uf 
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(29) 



(30) 



Thus for large U, Uf{U) — rif(U = oo) ~ U 2 . 
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FIGURE CAPTIONS 



Figure 1. (e +U/2)/\e \ vs. U for the symmetric Anderson lattice with V/2t=0. 375. e Q 
is the ground state energy per site. Our results are in good agreement with Gutzwiller 
and Monte Carlo |7j results. 



Figure 2a. Gaps vs. U for the symmetric case (t=l, V=l, e/ = —U /2). Open symbols 
are for the quasiparticle gap. Filled symbols are for the spin gap and the charge gap. 
Note that A s < A qp < A c for U > 0. 

Figure 2b. Gaps vs. U for the asymmetric case (t=l, V=l, e/=0). The charge and 
quasiparticle gap cross at U ~ 2. Note that the spin gap is smaller than both of them. 
Figure 3. Character of the spin and quasiparticle excitations vs. U . The open symbols 
are for < S* z > in the lowest excited S=l state at half filling. The filled symbols are for 
< Nf > —L in the ground state with 2L+1 electrons. The spin excitations have mainly 
f-character while the particle excitations have the character of conduction electrons in 
the strong coupling regime. 

Figure 4. Spin gap vs. 1/L 2 in the asymmetric case (t=l, V=l, e/=0, U=16). The 
intercept of the linear fitting gives the L = oo value. Inset: S z on a site vs. site for 
the lowest excited (S=l, S z = 1) state on a 24 site lattice at half filling. 
Figure 5a. f-spin-f-spin correlation function versus distance r apart for the symmetric 
case (t — 1, V — 1, 6f — —U/2, L=24). The oscillations increase in amplitude and 
persist for longer distances as U increases from 2 to 20. 

Figure 5b. f-spin-f-spin correlation function versus distance r apart for the asymmetric 
case (t — 1, V — 1, €f — 0, L — 24). As U increases from 2 to 20, neither the amplitude 
of the oscillations nor the correlation length change significantly. 
Figure 6. Correlation length vs. J e // = 8V 2 /U in the symmetric case. The curves for 
different sets of parameters fall on top of each other for J e ff < 1. The data for the 
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Kondo chain is taken from reference [36]. All data is for L=24. 

Figure 7. Staggered susceptibility vs. U for short chains with periodic boundary con- 
ditions (t=l, V=l). In the asymmetric case the susceptibility shows little change with 
U due to the suppression of the magnetic correlations in the ground state. Inset: Stag- 
gered susceptibility vs. q for short chains with U = 16. Notice that the asymmetric 
case (filled symbols) shows little dispersion indicating again that RKKY interactions 
are suppressed in the asymmetric case. 

Figure 8. The occupation rif of the f-level vs. U for L = oo in the asymmetric case 
(t=l, V=l, e-f = 0). As U — > oo, rif approaches a value close to 0.7. 
Figure 9a. log(A s ) vs. log(V) for the U = oo case (t=l, 6f = 0). The slave boson 
results agree qualitatively with our numerical results. 

Figure 9b. log(A s ) vs. log(V) for the U = oo case (t=l, e/ = —0.5). For small V, the 
slave boson results deviate significantly from our data. 

Figure 10a. log(l —nf) vs. log(V) for the U = oo case (t=l, ej = 0). The slave boson 
results agree qualitatively with our numerical results. 

Figure 10b. log(l — n/) vs. log(V) for the U = oo case (t=l, e/ = —0.5). For small 
V, the slave boson results deviate significantly from our data. 
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Figure 10a 
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Figure 10b 
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